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Abstract 

A lagrangian for relativistic fluid systems with matters inside is developed 
using gauge principle. In the model, the gauge boson represents the fluid 
field in a form = </>, where contains the fluid kinematics and (f> is 
an auxiliary field representing the fluid distribution. This leads to a new 
relativistic equation of motion for fluid, but which further coincides to the 
classical Euler equation at non-relativistic limit. The lagrangian is applied to 
model homogeneous universe as a bulk pure fluid system. Taking the simplest 
case of fluid with radial velocity and uniform distribution, the free energy 
density is calculated and its behaviour around Hubble distance is discussed. 

Keyword(s) : astrophysical fluid dynamics, gauge theory, lagrangian 
PACS : 98.90,+s; 98.80. Cq; ll.10.Ef; 11.15.-q 



*Email : lyman@tisda.org, sulaiman@teori.fisika.lipi. go. id 
^Email : tpdjun@teori.fisika.lipi. go. id 

*Email : handoko@teori.fisika.lipi. go. id, laksana.tri.handoko@lipi.go.id 

1 http://tcori.fisika.lipi.go.id 

2 http://tisda.bppt.go.id 

3 http:/ /www. fi.itb. acid 

4 http://www.fisika.ui.ac.id 



I Introduction 



Since long time ago, (especially relativistic) fluid dynamics has been applied widely 
in astrophysics to deal with some fluid like phenomena in our universe [1, 2, 3]. This 
approach is mainly motivated by the lack of interaction-based approach familiar 
in nuclear and particle physics that leads to in most cases unsolvable many body 
systems [4, 5]. The relativistic fluid is highly successful to describe the dynamics of 
many-body relativistic systems. It has been used to model some systems as small as 
heavy ion collision, and also as large as our universe with intermediate sized objects 
like neutron stars being considered along the way [6] . 

Most of those fluid-inspired models are constructed from the energy momentum 
tensor T^ v to take into account, not only the bulk motion of fluid, but also for ex- 
ample the random thermal velocity of each fluid particle and various forces between 
them which should contribute to the total potential energy of system. For instance, 
in the case of perfect fluid, it takes the well-known form T^ v = (p + P)-u M -u iy — Pg^ u , 
with density p, pressure P and = 7(1, v) for 7 = (1 — f3 2 ) ((3 = |v|/c). It is 
then argued that one can derive the continuity equation through the conservation 
of energy described by d^T^ = [7] . This method is actually analogous to classical 
fluid dynamics, that is starting from the continuity equation to derive the equation 
of motion (EOM) or vice versa. However, this kind of approach can not resolve 
the microscopic dynamics of a macroscopic bulk system modelled as a fluid system. 
This is crucial in some large scale objects like universe where inhomogeneity might 
significantly contribute to the whole system. 

In contrast with the classical approach, one can in principle construct the quan- 
tum fluid dynamics using the principle of quantum mechanics. For example, some 
works have been done in quantizing fluid either in the Hamiltonian framework using 
the Clebsch parametrization [8], or using canonical expressions [9, 10, 11]. Although 
the gauge symmetry has been adopted in these works, the models still require further 
generalization to extend them to the relativistic ones. There are also some works on 
gauging the fluid in the classical framework to incorporate electromagnetic interac- 
tion in a charged fluid system [12]. This kind of magnetohydrodynamics models and 
its relativistic versions are especially relevant for plasma physics as quark-gluon- 
plasma (QGP), heavy ion collision [13] and so on. 

Further advancements are adopting the lagrangian method. The main common 
advantage of this method is we are able to expand our theory to accomodate more 
physical interests like vorticity etc in a consistent manner since the initial properties 
of lagrangian would be conserved along the procedures [14, 15]. Moreover, the gauge 
extensions along these lines have been done for non-abelian fluid dynamics [16] , and 
for ideal fluid with global SO(3) symmetry [17]. We remark here that all of them 
have been developed starting from continuity equation. 

On the other hand, inspired by similarities between the Maxwell and fluid theo- 
ries [18], the classical fluid equation should be able to be reproduced directly from 
the electromagnetic lagrangian. However, in this paper we follow rather different 
approach, that is starting from a predefined gauge invariant relativistic lagrangian, 
and further deriving the EOM using Euler-Lagrange equation. This would avoid 
complexities on transforming the equations of non-relativistic to relativistic dynam- 
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ics, since we start from fully relativistic framework. Moreover, we propose a different 
form of gauge boson which might also be appropriate for performing numerical cal- 
culation using lattice gauge theory later on. 

The paper is organized as follows. First we introduce the model and derive the 
EOM and continuity equation for relativistic fluid. Before concluding the paper, we 
provide a simple example of modelling universe as a pure fluid system in the model, 
and calculate its free energy density. 



II The model 

Concerning the basic fluid properties which has no any intrinsic quantum numbers 
like spin, first let us consider a relativistic lagrangian (density) for bosonic matter, 

u = (d^y + v($) , (i) 

where V(<&) is the potential. For example in typical $ 4 theory, V(<&) = \m% + 

Inspired by the quantum electrodynamics (QED) theory, we put the lagrangian 
to be invariant under local Abelian gauge transformation [19, 20], U = exp[— i0{x)\ ~ 

1 — i6{x) with < 1. The boson field is then transformed as $ $' = 

exp[— id(x)\ $. It is well-known that the symmetry in Eq. (1) is revealed by in- 
troducing a gauge field which is transformed as —> Aj_ = A^ — ^(8^9), 
and replacing the derivative with the covariant one, = <9 M — ig A^. Further the 
gauge invariant kinetic term for gauge boson takes the form of F^ V F^ V , with strength 
tensor F^ u = d^A u - d u A^. 

Finally the total lagrangian with such gauge symmetry becomes, 

C = + C A + C mt , (2) 

where, 

Ca = -^V, (3) 
Ant = gJ^ + g 2 ($*$) A^A", (4) 
while the 4-vector current is, 

j^-it^-r^)] , (5) 

which satisfies the current conservation d^J^ = respectively. The coupling con- 
stant g represents the interaction strength between gauge field and matter. We note 
that the second term in Eq. (4) gives rise to the mass of fluid field, i.e. m 2 A = g 2 ($) 2 , 
where (<&) is the vacuum expectation value (VEV) of matter field. 

Rather than following the same argument in [18], here we construct the fluid 
lagrangian similar to the QED-like lagrangian as La above, that is considering A^ 
as a "fluid field" with velocity v. Further, let us take an ansatz that the gauge boson 
has the following form, 

A ll = {A , A) = (6) 
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with, 

e o = 7p 2 |v| 2 and a = -7^ , (7) 

where is an auxiliary boson field, while again 7 = (1 — f3 2 ) and (3 = |v|/c. We 
should remark that the auxiliary field is introduced in this paper to keep correct 
dimension. This actually represents the distribution of fluid in the system. Also, 
this separation is crucial in numerical calculation using lattice gauge theory later 
on, where we need to split the gauge field and its kinematic parameter to keep the 
velocity v appearing throughout the calculation. Because in lattice simulation the 
gauge boson is normally integrated out completely. However this point is out of 
coverage of this paper and will be discussed in detail in subsequent paper. 

Using Euler-Lagrange equation in term of against Eqs. (3) and (4), we obtain 
the following EOM, 

d,(d»A u ) - d 2 A, - gj, - g 2 A, = . (8) 

However, the last term in Eq. (8) breaks the gauge invariance. This can be resolved 
by taking the VEV of matter to be ($) = 0. So, the EOM becomes <9 M = gj^ 
where non-trivial relations are obtained for /i ^ v. We are now ready to derive the 
EOM relevant for relativistic fluid dynamics. Substituting Eqs. (6) and (7) into the 
EOM we obtain, 

for uniform 0, i.e. : constant. Here, Jj = J dxjJ = — f dtJj. Taking non- 
relativistic limit, i.e. (5 — > (7 — > 1), we find, 

^ + (v . V )v = -(vx W )- 9 j, (10) 

utilizing the vector identity |V |v| 2 = (v ■ V)v + v x (V x v), and uj = V x v 
is the vorticity. Actually, Eq. (10) reproduces the classical EOM for (turbulent) 
fluid. Therefore, one can conclude that Eq. (9) should be its relativistic version. 
Moreover, the lagrangian in Eq. (2) with A^ = e M written in Eq. (6) should 
describe a general relativistic fluid system interacting with the matters inside. 

We should remark few points here. First, the additional "current force" in Eqs. 
(9) and (10) is induced by the external source as manifestation of interacting fluid 
with the matters inside. Secondly, in the classical fluid dynamics the 4-vector cur- 
rent = (p, pv) represents the macroscopic distributions of fluid density (charge) 
and current. On the other hand, in our approach J M reveals the dynamics of the 
distribution function of matters inside the fluid as expressed in Eq. (5). This is 
analogous to the current in QED which is the result of fermion-pair interaction by 
emitting photon, and not the current of electromagnetic field itself as in the Maxwell 
theory. 

Now we are going to apply the lagrangian to model the universe. We should em- 
phasize that using the lagrangian would benefit us from calculating relevant physical 
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observables without dealing with unsolvable nonlinear EOM of fluid as shown in the 
subsequent section. Especially for gauge invariant lagrangian, we can further make 
use of known techniques in lattice gauge theory to perform any relevant calculations 
non-perturbatively. This would avoid unnecessary assumptions to enable perturba- 
tive calculations, since we should not in general assuming that g is small enough in 
the present model. 



Ill Example : pure and homogeneous fluid universe 



Let us consider about modelling a homogeneous universe without any matter inside 
as a relatisitic fluid system. This means the matter and interaction terms in the 
lagrangian can be omitted, and we take into account only the pure gauge term, that 
is C = La- 

In order to calculate, from the lagrangian, relevant physical observables in a 
system with finite temperature T in a unit volume, we bring the partition function 
density, 

Z = c-^ n , (11) 



where 7i is the hamiltonian density. Using T^ v = [dC/(d(d tJ- A p ))]d u A p — g pu C, 
we have the energy-momentum tensor being T pv = — (d p A p )F pu + \g pu F pa F pa . 
However, this is not symmetric and gauge invariant. The conserved, symmetric 



and gauge invariant T pu should be T pv 



-F pp F p u + \g pu F prj F pa by adding a term 
(d p A p )F pu . This additional term is admissible since in the case under consideration, 
no matter inside the fluid such that the current in Eq. (5) is absence. This leads 
to d p F pv = and then d p (A p )F pu = d p (A pj F pu ). Because the hamiltonian density is 
nothing else than the 00— element of T^, we obtain, 
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8 A „ , 



IV x A| 



(12) 



which is also already known in the electromagnetic theory as the electric and mag- 
netic terms. Assuming that the fluid has radial velocity, i.e. v = v(t,r) in polar 
coordinate, we calculate for each term, 



VxAl 



( dv dv\ 

= 0, 



(13) 
(14) 



by substituting A p in Eqs. (6) and (7), and putting uniform fluid distribution (0 : 
constant) in consequence of homogeneous fluid. 

The free energy density is given by T = —T \nZ. Therefore in the present case 
it yields, 



1 2 4 / dv 

r 7 U 



i0i 2 , 



(15) 



for steady velocity v(t,r) = v(r) using Eqs. (11), (12), (13) and (14). It is straight- 
forward to see that the free energy density is being infinite as (3 — > 1. In other 
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words, borrowing the Hubble law, v = rH [21], we can deduce from Eq. (15) that, 

T — > oo at r^r th = cH~ 1 . (16) 
This threshold distance r th is nothing else than the Hubble distance. 

IV Summary and discussion 

We have introduced a model for fluid universe with interacting matters inside. The 
lagrangian has been constructed using gauge principle and taking particular form 
of gauge boson as Eq. (6). The lagrangian should provide general description of 
relativistic fluid dynamics. This has been proven in the case of Abelian lagrangian 
with uniform field distribution which reproduces the classical Euler equation for 
turbulent fluid at non- relativistic limit. Using the lagrangian, we have calculated 
the free energy density for non-interacting homogeneous universe. In the present 
model the free energy would be getting infinite as approaching the Hubble distance, 
i.e. the universe horizon, where the system is completely relativistic. Beyond the 
horizon, the entropy is revived to be finite. 

According to this result, one might speculate that the horizon is like a mirror of 
twin universe if the fluid distributions in both sides are completely same. However 
one might take another functions for <fi with some reasons. For instance, if 4> follows 
a step function, 

*-«••)-{ \f ; :>:: ■ m 

with a constant a, the free energy beyond the horizon would be suppressed close to 
null for a quite large a. 

In contrary, at non-relativistic limit (/3 — > 0) the free energy should be propor- 
tional to the square of velocity. This result coincides with the classical kinetic theory 
of gas. 
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